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ABSTRACT 

We present the results of precision timing observations of the binary millisec- 
ond pulsar PSR J1640+2224. Combining the pulse arrival time measurements 
made with the Effelsberg 100-m radio telescope and the Arecibo 305-m radio 
telescope, we have extended the existing timing model of the pulsar to search 
for a presence of the effect of a general-relativistic Shapiro delay in the data. 
At the currently attainable precision level, the observed amplitude of the effect 
constrains the companion mass to m 2 = O.IS^qojMq, which is consistent with 
the estimates obtained from optical observations of the white dwarf companion 
and with the mass range predicted by theories of binary evolution. The measured 
shape of the Shapiro delay curve restricts the range of possible orbital inclinations 
of the PSR J1640+2224 system to 78° < i < 88°. The pulsar offers excellent 
prospects to significantly tighten these constraints in the near future. 

Subject headings: astrometry - stars: neutron - binaries: general - pulsars: 
individual (PSR J1640+2224) 
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1. 



INTRODUCTION 



Precision timing measurements of binary millisecond pulsars (Phinney & Kulkarni 1994) 
with sufficiently high (near edge-on) orbital inclinations make it possible to detect the effect of 
a general- relativistic time delay of the pulsar signal in the gravitational field of the companion 
star. For a pulsar in a circular orbit this "Shapiro delay" (Shapiro 1964) is given by 



where <£> is the orbital phase in radians, $o is the phase of the ascending node, and T = 
(G M Q /c 3 ). In practice, Shapiro delay is conveniently expressed in terms of two observables, 
the "range" r = m 2 T Q and the "shape" s = sini (Ryba & Taylor 1991), the post-Keplerian 
orbital parameters which allow a determination of the companion mass, m 2 , and the orbital 
inclination, i. 

Because high inclination orbits are relatively rare, Shapiro delay has been detected in 
only four pulsar-white dwarf (WD) binaries, PSR J1713+0747 (Camilo, Foster & Wolszczan 
1994), PSR B1855+09 (Ryba & Taylor 1991), PSR J0437-4715 (van Straten et al. 2001), 
PSR J1909-3744 (Jacoby et al. 2003), and possibly in PSR J0751+1807 (Nice, Splaver & 
Stairs 2003). Such detections offer very useful means to measure masses of the companion 
stars and to calibrate the pulsar spin-down models against the cooling models of white 
dwarfs. This is accomplished by comparing the spin-down age of a pulsar obtained from 
timing observations with the cooling age of a white dwarf estimated from its mass, an optical 
measurement of its temperature and from an appropriate cooling model (Kulkarni 1986). In 
particular, such comparisons are important in assessing the temperature modifying effect of 
hydrogen left over after the white dwarf formation (van Kerkwijk et al. 2000). 

Shapiro delay has also been detected in two double neutron star (NS-NS) systems, PSR 
B1534+12 (e.g. Stairs et al. 2002) and PSR J0737-3039A (Burgay et al. 2003; Lyne et al. 
2004). In these cases, the measured parameters r and s, together with the other two strong 
gravity related post-Keplerian parameters, the periastron advance u, and the time dilation 
and gravitational redshift 7, provide a "clean" test of general relativity and other theories of 
gravity, in the sense that it does not mix the relativistic strong-field and the radiative effects 
(Damour & Taylor 1992). 

Neutron stars in pulsar-WD binaries are thought to undergo extended periods of transfer 
of mass and angular momentum from their companions (Phinney & Kulkarni 1994). As a 
result, they are spun up to millisecond periods and may end up having masses significantly 
larger than the canonical value of 1.35 M derived for stars in the NS-NS systems (Thorsett 
& Chakrabarty 1999). The existing mass measurements for PSR J0437-4715 (mi = 1.58 ± 
0.18 M , van Straten et al. 2001) and PSR B1855+09 (m 2 = 1.57^ M ) (Nice, Splaver 
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& Stairs 2003, 2004) provide support for this idea. On the other hand, the recent mass 
measurement of 1.3 ±0.2 M derived from precision timing observations of PSR J1713+0747 
(Splaver et al. 2004) is in accord with the canonical value. Undoubtedly, more data for 
similar systems are needed to improve the existing statistics. 

In this paper, we report new results of timing measurements of the PSR J1640+2224 
binary system. A tentative detection of the Shapiro delay in the pulse arrival times from 
the pulsar allows us to set preliminary constraints on the orbital inclination of the system 
and the mass of the pulsar companion. In Sect. 2 we describe the timing observations 
made at Effelsberg and Arecibo, present the timing analysis, and summarize the resulting 
best-fit timing model for PSR J1640+2224. In Sect. 3 we discuss the new findings and 
their implications. In particular, we use simulated timing observations to demonstrate the 
expected potential of a future sub-microsecond timing of PSR J1640+2224 to verify the 
validity of the current best-fit model and to improve the estimates of masses of the pulsar 
and its white dwarf companion. 

2. OBSERVATIONS AND TIMING ANALYSIS 

We have conducted systematic, high-precision timing observations of PSR J1640+2224 
with the 100-m Effelsberg radiotelescope of the Max-Planck- Institut fur Radioastronomie 
in Bonn, Germany, and the 305-m Arecibo radiotelescope of the National Astronomy and 
Ionosphere Center in Puerto Rico, over a 7-year period from 1996 until 2003. 

At Effelsberg, PSR J1640+2224 was observed approximately once a month using a 
1300-1700 MHz tunable HEMT receiver at a centre frequency of 1410 MHz. In order to 
monitor changes of the dispersion measure (DM) we occasionally collected data at 860 MHz. 
As a backend, we used the Effelsberg-Berkeley Pulsar Processor (EBPP), which corrects for 
the dispersion smearing of the signal employing a coherent de-dispersion technique (Hankins 
& Rickett 1975). In the total power mode, the EBPP provided 32 channels for both senses 
of circular polarization with a maximum total bandwith of 112 MHz, depending on the DM 
and the observing frequency (Backer et al. 1997). For PSR J1640+2224, total bandwidths 
of 54 MHz and 27 MHz were available at 1410 MHz and 860 MHz, respectively. The output 
signals of each channel were fed into the de-disperser boards for coherent on-line de-dispersion 
and were synchronously folded at the pulse period over a 7 min integration time. 

At Arecibo, the timing observations of PSR J1640+2224 were made with the dual- 
circular polarization receiving systems at 430 MHz, 1130 MHz and 1410 MHz and the Penn 
State Pulsar Machine (PSPM). The PSPM pulsar backend is a computer controlled processor 



-4- 



with a 2 x 128 x 60 kHz filterbank designed to conduct fast sampled pulsar searches and 
precision timing measurements. Technical details of the backend are given in Cadwell (1997). 
For our timing observations, the two signals of opposite circular polarizations were added 
together, smoothed with a 32 /xs time constant, 4-bit quantized, folded synchronously with 
the topocentric pulse period, and stored for further processing. The pulse integration times 
were 3 min. at 430 MHz and 5 min. at both 1130 MHz and 1410 MHz. 

Both Effelsberg and Arecibo data were time stamped using the observatory hydrogen 
maser clocks and later synchronized to UTC(NIST) using the signals from the Global Po- 
sitioning System (GPS). In order to calculate the pulse time-of-arrival (TOA), high signal- 
to-noise template profiles of the pulse were constructed for each backend and observing fre- 
quency and least-squares fitted to the observed profiles in frequency domain (Taylor 1991). A 
theory- independent timing model for binary pulsars, devised by Damour & Deruelle (1986), 
was least squares fitted to the combined TOAs, weighted by their individual uncertainties, 
using the software package tempo 1 and the DE200 planetary ephemeris (Standish 1990). 

In the fitting procedure, the TOA segments obtained with the EBPP and the PSPM were 
fitted for an unknown offset between the two data sets resulting from different templates and 
TOA reference points in the profiles. Using the full TOA set at all frequencies we determined 
the DM of the pulsar. In the subsequent analysis, we fixed the best-fit value for the DM and 
used only the 1410 MHz TOAs from Effelsberg and both 1130 and 1410 MHz TOAs from 
Arecibo, as these high frequency data were not significantly affected by the observed DM 
variations. A determination of an initial best-fit model for the PSR J1640+2224 timing data 
involved a set of 12 parameters including the astrometric, and the rotational parameters of 
the pulsar and the orbital parameters of the binary system. In order to achieve a uniform 
reduced x 2 — 1 f° r eacn data segment, we increased the TOA uncertainties by a constant 
amount, approximately equal to the post-fit rms noise, by adding it in quadrature to the 
actual TOA values. 

In order to examine the timing data for a possible presence of the Shapiro delay we 
employed a grid search procedure used by Ryba & Taylor (1991). We searched the m 2 — 
cos? plane for a global x 2 minimum, by fixing the Shapiro parameters r and s at nodes of 
an appropriately defined 2D grid and repeatedly fitting for all other parameters for each 
set of (r, s) values. As displayed in Fig. 1, the grid search produces a well-defined global 
Ax 2 minimum, equivalent to the best-fit cos 2 = O.llioo? ~ &^-6 degrees) and vrt2 = 
0.151 

aos -^o (1°" uncertainties). Clearly, the inclusion of a Shapiro delay in the timing model 
for PSR J1640+2224 leads to astrophysically plausible estimates of both the companion 



1 http://pulsar.princeton.edu/tempo 
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mass and the inclination of the pulsar orbit. In Fig. 2, the timing residuals for the combined 
Arecibo and Effelsberg observations are plotted as a function of orbital phase. Because the 
observed Shapiro delay is weak and both r and s are strongly covariant with other model 
parameters, the effect is not detectable in residuals from the best fit involving the Keplerian 
orbit alone, as seen in Fig. 2a. On the other hand, in Fig. 2b, showing a Shapiro delay 
signature extracted from the grid search with all other timing effects removed (see also Ryba 
& Taylor (1991) and Camilo, Foster & Wolszczan (1994)), the amplitude of the effect in the 
PSR J 1640+2224 TOA residuals significantly exceeds the TOA uncertainties, as expected 
from the result of the \ 2 search displayed in Fig. I. 

In principle, the observed signature could be induced by DM variations over the pul- 
sar orbit. If the pulsar's WD companion had an extended envelope created by the pulsar 
wind and the high-energy photon flux, as observed in eclipsing binary systems (e.g. Nice, 
Arzoumanian & Thorsett 2000), the electron column density would fluctuate periodically as 
a function of orbital phase. For highly inclined orbits this would obviously cause periodic, 
frequency-dependent TOA variations that could mimic the effect of Shapiro delay. We have 
ruled out this possibility by verifying that the effect has the same amplitude in the TOA 
measurements made at four different frequencies. 

In the case of binary pulsars with nearly circular, low-inclination orbits, Shapiro delay 
becomes covariant with Roemer delay and cannot be measured (Lange et al. 2001). PSR 
J1640+2224 has the most eccentric orbit among the pulsar-WD binaries (see Table 1 and 
Edwards and Bailes 2001) and its inclination angle of i = 84 ± 6° derived from our analysis 
appears to be high enough to allow the inclusion of Shapiro delay in the timing model 
(Fig. 2b). In any case, further observations of the pulsar with higher timing precision are 
clearly necessary to fully assess a statistical significance of our detection. 

The parameters of the best-fit timing model for PSR J1640+2224 are listed in Table 1 
along with the ones for which only upper limits could be determined. In this case, the upper 
limits were obtained by allowing the parameters to vary, one at a time, in the global fit. 
Also included in the table are the most important parameters derived from the final model. 
Finally, the behavior of the post-fit timing residuals as a function of time, spanning a 7-year 
period, is shown in Fig. 3. Evidently, the best-fit model including the Shapiro delay leaves 
no additional systematic effects above the current post-fit rms residual of 2.0 /is. 
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3. DISCUSSION 

A new timing model for the binary millisecond pulsar PSR J1640+2224 discussed in this 
paper is entirely consistent with the previous models published by Wolszczan et al. (2000) 
and Potapov et al. (2003). In addition, owing to a higher timing precision, the new model 
provides further reduction of the parameter estimation errors and, above all, it includes 
astrophysically sensible estimates of the Shapiro delay parameters. 

The binary companion to PSR J1640+2224 is a white dwarf with the estimated cooling 
age and mass of 7±2 Gyr and 0.25±0.10M Q , respectively, as determined from the Hubble 
Space Telescope observations (Lundgren et al. 1996). A range of masses predicted by the 
binary period-companion mass {P b — m 2 ) relationship based on the theory of low- and 
intermediate-mass binary evolution is 0.35 < m 2 < 0.39 M (Tauris & Savonije 1999), and 
the minimum companion mass from the mass function, /(mi, 777.2) = (^2 sini) 3 (mi+m 2 ) _2 = 
(2n/P b ) 2 x 3 /T Q = 0.0059 M Q is m 2 = 0.25 M for a m x = 1.4 M neutron star. At the 
presently attainable level of accuracy, the best-fit companion mass of m 2 = 0.151q 05-^0 
derived from our data is consistent with the above estimates, as illustrated in Fig. 1. 

Among the pulsar- WD binaries with a detectable Shapiro delay, only PSR J0437— 4715, 
PSR J1713+0747, and PSR B1855+09 have the values of m 2 and sin 2 determined with the 
accuracy that is high enough to make them usable in setting a tight constraint on the pulsar 
mass (van Straten et al. 2001; Nice, Splaver & Stairs 2003, 2004; Splaver et al. 2004) and 
in investigating the details of evolution of the pulsar's WD companion (van Kerkwijk et 
al. 2000). We have examined a future potential of the PSR J1640+2224 timing to become 
comparably useful by generating artificial TOAs according to the model of Table 1 and 
analyzing the data over progressively longer periods of time for several reasonable values of 
the timing precision. Encouragingly, as shown in Fig. 4, in only four years of monthly timing 
measurements with a 0.5 fis precision, the estimation errors of m 2 and cos 2 approach the 
respective levels of 0.01 M Q and 0.001. Since the Arecibo timing measurements using the 
PSPM and an 8 MHz receiver bandwidth are characterized by a ~1 /is long-term residual, 
it is quite conceivable that the required <0.5 /is precision can be achieved with a new 
generation of broadband, 100 MHz bandwidth backends already available at the telescope. 
Further observations at this level of precision will quite conceivably allow a verification of the 
timing model presented in this paper. PSR J1640+2224 is also likely to become a valuable 
member of a set of the most accurate pulsar clocks that can be timed either individually, or 
as an array to detect a low-frequency background of gravitational radiation (e.g. Thorsett 
k Dewey 1996, Jaffe & Backer 2003). 

We are very grateful to all staff at the Effelsberg and Arecibo observatories for their 
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i , 1 , 1 , 1 1 

0.2 0.4 0.6 

cos i 

Fig. 1. — Constraints on m 2 and cost in the PSR J1640+2224 system from a x 2 search for 
the best-fit Shapiro delay parameters. The global x 2 minimum is indicated by a cross. The 
contours of A% 2 = 1.0 (solid), A% 2 = 4.0 (dashed-dotted) and A% 2 = 9.0 (dashed) have 
extrema respectively corresponding to la, 2a, and 3a errors on the individual parameters 
cosi and m 2 . Lines of constant mi are indicated. Horizontal lines denote the m 2 -range 
bounded by a lower limit obtained from optical observations (m 2 > 0.15M Q ) and an upper 
limit allowed by the P b — m 2 relationship (m 2 < 0.39M Q ). 
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Fig. 2. — Timing residuals for PSR J1640+2224 observed with the Arecibo telescope at 1130 
MHz and 1410 MHz (filled circles) and with the Effelsberg telescope at 1410 MHz (open 
circles), as a function of orbital phase, (a) Post-fit residuals for the best-fit model involving 
only the five Keplerian parameters, (b) The effect of Shapiro delay on the TOA residuals 
calculated with sku and m 2 set to zero and all other parameters fixed at their best-fit values. 
The solid curve represents the delay predicted by general relativity for the best-fit Shapiro 
parameters, (c) Post-fit residuals for the best-fit model including the sin i and m 2 parameters 
(see also Table 1). 
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Fig. 3. — Best-fit timing residuals for PSR J1640+2224 as a function of observing epoch. 
See Fig. 2 for further explanation. 
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Fig. 4. — Estimation errors of the companion mass and orbital inclination from the simulated 
TOA measurements of PSR J1640+2224 with a 0.5 /is precision. The initial error values are 
equal to those currently observed. 
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Table 1: Timing model for PSR J1640+2224 



Measured parameters 
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Epoch (MJD) 


51700.0 


Dispersion measure, DM (pccm -3 ) 


18.4260(8) 


Orbital period, Pb (d) 


175.46066194(7) 


Projected semi-major axis, x (lt-s) 


55.3297198(4) 


Eccentricity, e 


0.000797262(14) 


Epoch of periastron 6 , T (MJD) 


51626.1785(5) 


Longitude of periastron b , u (deg) 


50.7308(10) 


Shape of Shapiro delay c , s 


n qq+°- 01 


Range of Shapiro delay c , r (fj,s) 


n 74+0.39 

u - ' ^-0.25 


Measured upper limits' 1 


Parallax, n (mas) 


< 3.7 


Pulse frequency second derivative, (s~ 3 ) 


< 4 x IO" 27 


DM derivative, \d(DM)/dt\ (pccm^yr- 1 ) 


< 1.3 x 10~ 3 


Orbital period derivative, Pb (ss -1 ) 


< 3 x IO" 10 


Derivative of projected semi-major axis, \x\ (lt-s s^ 1 ) 


< 1.7 x 10~ 14 


Periastron rate of change, \u\ (deg yr' 1 ) 


< 1.1 x IO" 3 


Derived parameters 


Galactic longitude, / 


41?051 


Galactic latitude, b 


38?271 


DM distance 6 (kpc) 


1.16 


Composite proper motion, \i (masyr -1 ) 


11.4(2) 


Companion mass, m-i (M ) 


o.i5l8:8i 


Orbital inclination angle, i (deg) 


84+ 4 


Mass function, / m (M Q ) 


0.0059074304(2) 


Number of TOAs 


314 


Timing RMS (jib) 


2.0 



a Figures in parentheses are 2a uncertainties in the last digits quoted (twice the formal tempo 
errors). 

b lo and To are highly covariant. Observers should use u> = 50.730834835740 and 
T = 51626.178534099. 

c Uncertainties are la errors derived from the \ 2 analysis (see Sect. 2). 



